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Abstract 



We present the results of a search for anomalous production of dipho- 
ton events with large missing transverse energy using the Collider Detec- 
tor at Fermilab. In 202 pb _1 of pp collisions at \fs = 1.96 TeV we ob- 
serve no candidate events, with an expected standard model background of 
0.27 ± 0.07(stat) ± O.lO(syst) events. The results exclude a lightest chargino 
of mass less than 167 GeV/c 2 , and lightest neutralino of mass less than 
93 GeV/c 2 at 95% C.L. in a gauge-mediated supersymmetry-breaking model 
with a light gravitino. 



PACS numbers 13.85Rm, 13.85Qk, 14.80.-j,14.80.Ly 
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The standard model (SM) [1] of elementary particles has been enormously successful, 
but it is incomplete. For theoretical reasons [2,3], and because of the 'ee77+missing trans- 
verse energy (faY [4] candidate event recorded by the CDF detector in Run I [5], there is 
a compelling rationale to search in high-energy collisions for the production of heavy new 
particles that decay producing the signature of 77 + fa. Of particular theoretical interest 
are supersymmetric (SUSY) models with gauge-mediated SUSY-breaking (GMSB). Char- 
acteristically, the effective SUSY-breaking scale (A) can be as low as 100 TeV, the lightest 
SUSY particle is a light gravitino (G) that is assumed to be stable, and the SUSY particles 
have masses that may make them accessible at Tevatron energies [2]. In these models the 
visible signatures are determined by the properties of the next-to-lightest SUSY particle 
(NLSP) that may be, for example, a slepton or the lightest neutralino (xi)- I n the GMSB 
model investigated here, the NLSP is a Xi decaying almost exclusively to a photon and a G 
that penetrates the detector without interacting, producing fa. SUSY particle production 
at the Tevatron is predicted to be dominated by pairs of the lightest chargino (xf) an d by 
associated production of a xt an d the next-to-lightest neutralino (X2). Each gaugino pair 
cascades down to two x?'s, leading to a final state of 77 + fa + X, where X represents any 
other final state particles. 

In this paper we summarize [6] a search for anomalous production of inclusive 77+^+X 
events in data corresponding to an integrated luminosity of 202 ± 12 pb _1 [7] of pp collisions 
at y/s = 1.96 TeV using the CDF II detector [8]. We examine events with two isolated 
photons with \rj\ < 1.0 and E T > 13 GeV for the presence of large fa. This work extends a 
previous CDF search [5] for SUSY in this channel by using an upgraded detector, a higher 
pp center-of-mass energy, and a larger data sample. The analysis selection criteria have been 
re-optimized to maximize, a priori, the expected sensitivity to GMSB SUSY based only on 
the background expectations and the predictions of the model. Similar searches for diphoton 
+ fa events have been performed elsewhere [9]. 

We briefly describe the aspects of the CDF II detector relevant to this analysis. The 
magnetic spectrometer consists of tracking devices inside the 3-m diameter, 5-m long su- 
perconducting solenoid magnet operating at 1.4 T. A 90-cm long silicon micro-strip vertex 
detector, consisting of one single-sided layer and six double-sided layers, with an additional 
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double-sided layer at large 77, surrounds the beam pipe. Outside the silicon detector, a 3.1-m 
long drift chamber with 96 layers of sense wires is used with the silicon detector to deter- 
mine the momenta of charged particles and the z position of the pp interaction (2 vertex ). The 
calorimeter, constructed of projective towers, each with an electromagnetic and hadronic 
compartment, is divided into a central barrel that surrounds the solenoid coil < 1.1) 
and a pair of 'end-plugs' that cover the region 1.1 < \r)\ < 3.6. The hadronic compartments 
of the calorimeter are also used to provide a measurement of the arrival time of the par- 
ticles depositing energy in each tower. Wire chambers with cathode-strip readout (the 
CES system), located at shower maximum in the central electromagnetic calorimeter, give 
2-dimensional profiles of showers. A system of proportional wire chambers in front of the 
central electromagnetic calorimeters (the CPR system) uses the one-radiation-length-thick 
magnet coil as a 'preradiator' to determine whether showers start before the calorimeter [10]. 
Muons are identified with a system of planar drift chambers situated outside the calorimeters 
in the region \r]\ < 1.0. 

We select candidate events using both online (during data taking) and offline selection re- 
quirements. Online, events are selected for the presence of two photon candidates, identified 
by the three-level trigger as two isolated electromagnetic clusters [10] with E% > 12 GeV, 
or two electromagnetic clusters with Ej. > 18 GeV and no isolation requirement. The of- 
fline event selection requirements for the diphoton candidate sample are designed to reduce 
electron and jet/7r° backgrounds while accepting well- measured diphoton candidates. We 
require two central (approximately 0.05 < \r)\ < 1.0) electromagnetic clusters that: a) have 
Ej, > 13 GeV; b) are not near the boundary in of a calorimeter tower [11]; c) have the ratio 
of hadronic to electromagnetic energy, Had/EM, < 0.055 + 0.00045 • -E 7 (GeV _1 ); d) have 
no tracks, or only one track with < 1 GeV/c, extrapolating to the towers of the cluster; 
e) are isolated in the calorimeter and tracking chamber [12]; f) have a shower shape in the 
CES consistent with a single photon; g) have no other significant energy deposited nearby 
in the CES. 

To minimize the number of events with large fir due to calorimeter energy mis- 
measurement, we correct for jet (J) energy loss in cracks between detector components 
and for nonlinear calorimeter response [13]. To avoid any remaining cases where a jet is not 
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fully measured by the calorimeter, we remove events based on the azimuthal opening angle 
between the fir direction and the of any jet with uncorrected E T > 10 GeV, A<j)(fir,j). 
We require all events to have 10° < A<j)(fir, j) < 170°. To reduce beam-related and cosmic- 
ray backgrounds we require a good vertex with |-2 vertC x| < 60 cm and reject events with 
significant energy out-of-time with the collision [14]. These backgrounds can also produce 
fir equal in magnitude and opposite in direction to a photon, or to the vector sum of the 
momenta of two photons if they are nearby in 0. In this case an event is rejected if there 
are potential cosmic-ray hits in the muon chamber, within 30 degrees of the photon, that 
are not matched to any track. Events are also rejected if there is a pattern of energy in 
the calorimeter indicative of beam-related backgrounds [15]. A sample of 3,306 diphoton 
events pass all candidate selection requirements. The fir requirement, fir > 45 GeV, is de- 
termined by the final optimization procedure that is discussed below, after a more complete 
description of the backgrounds. 

Before the fir requirement, the diphoton candidate sample is dominated by QCD inter- 
actions producing combinations of photons and jets faking photons. In each case only small 
measured fir is expected, due mostly to energy measurement resolution effects. Standard 
CDF techniques [10] are used to estimate the individual contributions for the sample to 
be 47 ± 6% 7j, 29 ± 4% 77, and 24 ± 4% jj production. To estimate the shape of the 
fir distribution of this background we use a control sample of similarly-produced events 
that have the same calorimetric response and resolution. We select 7,806 events that pass 
the same photon E T , z VCTtcx , fiducial, A(fi(fir,j), beam-related and cosmic-ray background 
selection requirements, but are allowed to satisfy looser photon identification and isolation 
requirements [16]. If an event is in the diphoton candidate sample it is rejected from the 
control sample. The contribution from events, discussed below, is also subtracted from 
the control sample. Since the fir resolution for a given event is a function of the sum of 
all the transverse energy in the event (Y,E T ), and we observe a small difference between 
the T,E T distributions of the diphoton candidate and control samples, we correct the fir in 
the control sample for this difference [17]. To predict the number of events with large fir, 
we normalize the corrected control sample distribution to the number of diphoton candi- 
date events in the region fir < 20 GeV, and fit the spectrum above 10 GeV to a double 
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exponential. We predict 0.01 ± O.Ol(stat) ± O.Ol(syst) events with fir > 45 GeV, where 
the uncertainty is dominated by differences in the predictions using various control sample 
selection requirements, the choice of fit function, and the statistical uncertainties of the 
sample. 

Events with an electron and a photon candidate (W7 — > ez/7, Wj — > eu^f a k ei -Z7 — ■> ee7, 
etc.) can contribute to the diphoton candidate sample when the electron track is lost (by 
tracking inefficiency or bremsstrahlung) to create a fake photon. For W decays large fir can 
come from the neutrinos. This background is estimated using events from the data. The 
diphoton triggers accept electromagnetic clusters with tracks so they provide an efficient 
and unbiased sample of these events. We find 462 &y events before the fir requirement. 
Examining a Z — > ee sample, we estimate 1.0 ± 0.4% of electrons will pass the diphoton 
candidate sample requirements, including charged track rejection. By multiplying the num- 
ber of observed e^fir events by the probability that an electron fakes a photon, we estimate 
0.14 ± 0.06(stat) ± 0.05(syst) background events in the sample with fir > 45 GeV. The 
uncertainty is dominated by the statistical uncertainty in the fake rate and the uncertainty 
in the purity of the sample. 

Beam-related sources and cosmic rays overlapped with a SM event can contribute to 
the background by producing spurious energy deposits that in turn affects the measured 
fir- While the rate at which these events contribute to the diphoton candidate sample is 
low, most contain large fir- The spurious clusters can pass photon cuts. The dominant 
contribution actually comes from sources that produce two photon candidates at once, such 
as a cosmic muon undergoing bremsstrahlung twice. This background is estimated from the 
data using a sample of events with no primary collision and two electromagnetic clusters, 
multiplied by the rate that clusters from cosmic rays pass the diphoton candidate sample 
requirements. Backgrounds where only one of the photons, or only the fir, is from a non- 
collision source, are estimated to be negligible. The total number of events expected from 
non-collision sources in the fir > 45 GeV sample is 0.12 ± 0.03(stat) ± 0.09(syst). The 
uncertainty includes the uncertainty in the rate that spurious clusters pass the diphoton 
selection requirements and takes into account the statistics and purity of the sample of 
events with no primary collision. 
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The pr distribution of the diphoton candidate sample, see Figure 1, shows good agree- 
ment with that from the expected backgrounds. Table I summarizes the number of observed 
events and predicted backgrounds with four different fir requirements. There are no events 
with fir > 45 GeV. 

Since there is no evidence for events with anomalous fir in the diphoton candidate 
sample, we set limits on new particle production from GMSB using the parameters suggested 
in Ref. [18]. To estimate the acceptance for this scenario we generate GMSB events using 
ISAJET [19] with CTEQ5L parton distribution functions [20]. The production cross sections 
from ISAJET are corrected by a .fT-factor of approximately 1.2 to match the next-to- leading 
order (NLO) prediction [21]. We process the events through the GEANT-based [22] detector 
simulation, and correct the resulting efficiency with information from data measurements. 

Since electrons and photons interact similarly in the calorimeter we investigate the effi- 
ciency of the photon identification and isolation selection criteria by using a control sample 
of electrons from Z — » ee events. Separate efficiency estimates comparing data and detector 
simulation agree to within 3%. Using the simulation we estimate that if a photon within 
the fiducial portion of the detector is isolated, it has an 80% probability of passing the iden- 
tification and isolation criteria. However, the isolation energy of the photons is predicted 
from the Monte Carlo to be a strong function of the SUSY scale due to the number and 
energy of the extra jets produced. We find, for example, the single-photon efficiency to 
be reduced to 62% at M~± = 170 GeV/c 2 . This has a significant impact on the sensitivity. 
We find that the fraction of generated signal events passing all the selection requirements, 
including fir > 45 GeV, rises linearly from 3.5% at M~±=100 GeV/c 2 to approximately 8% 

Al 

at 180 GeV/c 2 . It remains roughly flat for larger masses due to the increasing inefficiency 
of the A<p(fir,j) selection requirement. The relative systematic uncertainty in the efficiency 
of the photon identification and isolation requirements is approximately 6.5% per photon. 
Other significant uncertainties in the Monte Carlo model predictions are from initial/final 
state radiation (10%), Q 2 of the interaction (3%) and uncertainty in parton distribution 
functions (5%). Combining these numbers with the 6% luminosity uncertainty gives a total 
relative systematic uncertainty of 18%. 

The kinematic selection requirements defining the final data sample are determined by 
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a study to optimize the expected limit, i.e., without looking at the signal region data. To 

compute the expected 95% confidence level (C.L.) cross section upper limit we combine 

the predicted signal and background estimates with the systematic uncertainties using a 

Bayesian method [23] and follow the prescription described in Ref . [24] . The expected limits 

are computed as a function of fir, photon Et, and A<f>(fir,j) selection requirements. We find 

that the best limit is predicted with the selection described above for the diphoton candidate 

sample, and fir> 45 GeV. The statistical analysis indicates that the most probable expected 

result, in the absence of a signal, would be an exclusion of M~± less than 161 GeV/c 2 and 

Xi 

M~o less than 86 GeV/c 2 . 
Xi 

In the data signal region, with fip> 45 GeV, we observe zero events. Taking into account 
the 18% systematic uncertainty we set a 95% C.L. upper limit of 3.3 signal events. Figure 2 

shows the observed cross section limits as a function of M~± and M~o along with the 

Xi Xi 

theoretical LO and NLO production cross sections. Using the NLO predictions we set a 

limit of M~± > 167 GeV/c 2 at 95% C.L. From mass relations in the model, we equivalently 

Xi 

exclude M~ < 93 GeV/c 2 and A < 69 TeV. 
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FIGURES 
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FIG. 1. The fJr spectrum for events with two isolated central photons with Ej, > 13 GeV 
and \r]\ < 1.0 along with the predictions from the GMSB model with a xf mass of 175 GeV/c 2 , 
normalized to 202 pb _1 . The diphoton candidate sample data are in good agreement with the 
background predictions. There are no events above the pr > 45 GeV threshold. The properties of 
the two candidates above 40 GeV appear consistent with the expected backgrounds. 
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FIG. 2. The 95% C.L. upper limits on the total production cross section times branching ratio 

versus M~± and M~o for the light gravitino scenario using the parameters proposed in [18]. The 
Xi Xi ' 

lines show the experimental limit and the LO and NLO theoretically predicted cross sections. We 
set limits of M~± > 167 GeV/c 2 and M~ > 93 GeV/c 2 at 95% C.L. 

Xl Al 

In conclusion, we have searched 202 pb _1 of inclusive diphoton events at CDF run II 
for anomalous production of missing transverse energy as evidence of new physics. We find 
good agreement with standard model expectations. We find no events above the a priori 
frp threshold, and thus observe no new ee^jftr candidates. Using these results, we have 
set limits on the lightest chargino M~± > 167 GeV/c 2 and M~ > 93 GeV/c 2 at 95% C.L. 

Al Al 

in a GMSB model. This limit is an improvement over previous CDF and D0 limits and is 



13 



comparable to LEP II for similar models [9]. 
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TABLES 





Expected 


Observed 


Requirement 


QCD 




Non-Collision 


Total 




25 GeV 


4.01 ±3.21 ±3.76 


1.40 ±0.52 ±0.45 


0.54 ±0.06 ±0.42 


5.95 ±3.25 ±3.81 


3 


35 GeV 


0.30 ±0.24 ±0.22 


0.84 ±0.32 ±0.27 


0.25 ±0.04 ±0.19 


1.39 ±0.40 ±0.40 


2 


45 GeV 


0.01 ±0.01 ±0.01 


0.14 ±0.06 ±0.05 


0.12 ±0.03 ±0.09 


0.27 ±0.07 ±0.10 





55 GeV 


(negligible) 


0.05 ±0.03 ±0.02 


0.07 ±0.02 ±0.05 


0.12 ±0.04 ±0.05 






TABLE I. Numbers of events observed and events expected from background sources as a 
function of the pr requirement. Here "QCD" includes the 77, 7J and jj processes. The first 
uncertainty is statistical, the second is systematic. 
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